Complexity and line of critical points in a short-range spin-glass model 
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We investigate the critical behavior of a three-dimensional short-range spin glass model in the 
presence of an external field e conjugated to the Edwards- Anderson order parameter. In the mean- 
field approximation this model is described by the Adam-Gibbs-DiMarzio approach for the glass 
transition. By Monte Carlo numerical simulations we find indications for the existence of a line 
of critical points in the plane (e, T) which separates two paramagnetic phases and terminates in a 
critical endpoint. This line of critical points appears due to the large degeneracy of metastable states 
present in the system (configurational entropy) and is reminiscent of the first-order phase transition 
present in the mean-field limit. We propose a scenario for the spin-glass transition at e = 0, driven by 
a spinodal point present above T c , which induces strong metastability through Griffiths singularities 
effects and induces the absence of a two-step shape relaxation curve characteristic of glasses. 



Among all different approaches to understand the glass 
transition the thermodynamic theory of Adam-Gibbs-Di 
Marzio (AGM) has deserved a lot of interest during the 
last decades 0]. The AGM theory predicts the occur- 
rence, at a Kauzmann temperature Tk, of a second or- 
der phase transition for the metastable undercooled liq- 
uid where the configurational entropy (also called com- 
plexity) vanishes. The validity of the AGM theory for 
real glasses has never been demonstrated so the correct 
description of the glass transition still remains open [^J. 
An alternative dynamical approach was proposed in the 
eighties to describe relaxational processes in the under- 
cooled liquid regime, experimentally observed in scatter- 
ing and dielectric measurements. This received the name 
of mode-coupling theory (MCT) ||. 

Quite recently it has been realized that spin glasses 
are models which account for both the thermodynamic 
(AGM) and the dynamical (MCT) approaches j|. Al- 
though spin glasses are models with quenched disorder 
(and structural glasses are not) this is not an essential 
difference because the existence of a crystal phase in 
structural glasses has no dramatic effect in the dynamics 
of the (disordered) metastable undercooled liquid phase. 
Unfortunately, up to now this connection between spin- 
glasses and glasses remains only at the mean-field level 
and it is not clear what happens if one considers short- 
range interactions. In fact, concepts such as complexity 
in AGM or the ergodicity parameter in ideal MCT are 
originally mean-field and it is not clear what is their rel- 
evance in short-ranged realistic systems. 

Recently it has been suggested |||| that the effect of 
the complexity could be observed through numerical sim- 



ulations in a generic glassy system coupling two replicas 
by introducing a term in the Hamiltonian of the type 
—eq (e being the conjugate field of the order parameter 
q which is the overlap between the configurations of the 
two replicas ) . Through the study of an exactly solv- 
able spin-glass model it has been shown the existence of 
a first order transition line T c (e) with a critical end-point 
Q . Their result is an explicit check of the fact that the 
glass transition for e — (where the complexity van- 
ishes) is a first order phase-transition (in the sense that 
the order parameter q is discontinuous) and the point 
T c = T c (e = 0) is a tricritical point. Again, this result 
has been obtained within the mean-field approximation 
and it is unclear to what extent this result is valid in a 
finite-dimensional model. Recent numerical simulations 
on a short-range version of p-spin Ising spin glass 
have shown that the mean-field discontinuous transition 
becomes continuous in finite dimensions . So the charac- 
teristic first-order transition predicted within mean-field 
theory dramatically changes in finite dimensions. In this 
work we want to show that this is not completely true 
and that some features of the mean-field approximation 
indeed survive in finite dimensions although they must 
be appropriately interpreted. This allows to interpret our 
findings in terms of a picture for a continuous spin-glass 
transition induced by the strong metastability and driven 
by the collapse of the complexity or configurational en- 
tropy. Furthermore, our results point in the direction 
that disordered systems in three dimensions are very well 
described by a line of critical points (characteristic of sys- 
tems at their lower critical dimension) in agreement with 
recent numerical simulations in the Edwards-Anderson 
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model in three dimensions |10|. 
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FIG. 1. Different phases for the Edwards- Anderson model 
(above) and the short-range p-spin model (below). In the for- 
mer case there is a single paramagnetic phase, in the latter 
two different paramagnetic phases divided by a critical line 
which connects a tricritical point (TCP) and a critical end- 
point (CE) (for p even the line TP — CE exists also for e < 0). 
In both cases the spin-glass phase is present only for e = 0. 

We shall consider a recently introduced short-range p- 
spin glass model |8]|J which is defined on a d dimensional 
hypercubic lattice. On each site of the lattice there are 
M spins interacting with the following Hamiltonian 
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By V^- .we sum over all sites of the lattice all 
possible groups of p spins that can be formed between 
spins on the same site and spins of adjacent sites. In this 
work we consider fsing spin variables. By <Ji r lr we denote 
the l^ 1 spin of site i r with the index l r running from 1 to 
M . For p — 2 ([!]) corresponds to the Edwards- Anderson 
model. Although this model may seem quite artificial it 
has the advantage of including multispin interactions in 
a lattice without introducing new symmetries which may 
change the degeneracy of the ground state | Q . 

Here we consider two identical coupled models (each 
defined through (|l|)) via the following Hamiltonian 



n 2 p (M, M) = n p ({a}) + n p ({r}) - ev q 
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where V = L 3 is the volume and Vq 
fines the order parameter. Because the e = transition in 
this model is continuous we naively expect that (similarly 



as it happens in the Edwards- Anderson model) there is 
no phase transition for e ^ 0. In mean-field models with 
a one-step replica broken phase, a transition exists for 
finite e 0. In that case, there is a transition line which 
separates two paramagnetic phases with a finite latent 
heat (which vanishes at the critical endpoint). The two 
possible phases are depicted in figure 0. The Edwards- 
Anderson order parameter qEA (defined as the infinite- 
time limit of the equilibrium autocorrelation function H] ) 
displays a finite jump across the line T c (e) which vanishes 
at the critical endpoint. Here we find strong indications, 
through Monte Carlo simulations, that this first-order 
line T c (e) persists in finite dimensions but becomes a line 
of critical points. So, in finite dimensions the first-order 
line becomes continuous (i.e. qEA is continuous when 
crossing the transition line and there is no latent heat), 
the critical endpoint displaying a higher-order singular- 
ity. 

In this paper we will focus our research of ([!]) for case 
M = 3,p = 4 in D = 3. Measurements of the spin-glass 
susceptibility for ([!]) show that this model has a continu- 
ous finite-temperature phase transition at T c ~ 2.6 with 
a divergent spin-glass susceptibility and a small nega- 
tive specific heat exponent ||. To evidenciate a phase 
transition for finite e which separates two different para- 
magnetic phases we have done a detailed Monte Carlo 
study |T^] of the Binder parameter as a function of both 
the coupling e and the temperature in the paramagnetic 
region. The Binder parameter is usually defined through 
the relation, 



ff(e,T) = i(3- 



(g-<g>) 4 ■ 

2\2 - 



((q-<q>) 2 ) 



(3) 



where < .. > stands for statistical average and (.) for 
disorder average. Concerning the Binder parameter we 
expect that it should vanish everywhere in both param- 
agnetic phases except at the critical line where it should 
be finite. So, if we fix e (or, equivalently the tempera- 
ture) and vary the temperature (equivalently e) we expect 
the presence a maximum approximately located on the 
transition line at a temperature T C (L, e). The results for 
g(e,T) are shown in figure || for T — 3.2 and different 
sizes. Note that the Binder parameter shows a maxi- 
mum located at e = 0.05. As a comparison we have done 
simulations for the Edwards- Anderson model in three di- 
mensions above T c by coupling two replicas which evi- 
denciate the absence of a maximum in this case. So the 
presence of a maximum in g(e,T) already for small sizes 
is a main feature of this model. For L = 3, 4, 5 finite-size 
corrections appear to be quite strong (this was already 
observed in Q by measuring the P{q)) and the position 
of the maximum of the Binder parameter as well as its 
value both shift with L. Nevertheless, the maximum of 
the Binder parameter for L = 6 superimposes with the 
maximum for L — 5 and the Binder parameter goes to 
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zero outside from the maximum as L increases. This 
result indicates the presence of a critical line which sep- 
arates two paramagnetic phases. A more stringent test 
of this result requires simulating larger sizes than those 
we did. Unfortunately for L = 7 the thermalization time 
is much larger than what we can afford with the present 
numerical methods. Actually, a test of the needed ther- 
malization time shows that it grows dramatically with 
L and for L = 7 and T — 3.2 this may be larger than 
several hundreds of million of Msteps for a non vanishing 
fraction of disorder realizations. As we will argue later 
this is consequence of the strong metastability and the 
highly corrugated landscape characteristic of this model 
which, on the other hand, induces the existence of this 
critical line (actually thermalization is much easier in the 
Edwards- Anderson model). Note that the fact that the 
maximum value of g along the critical line is smaller than 
one (approximately 0.2) is an indication that the transi- 
tion is continuous in q (if there were a finite jump in q the 
maximum value of g should be 1). Moreover, our data 
do not show any indication of a jump in the value q(e) 
as a function of the temperature (or q(T) as a function 
of e) in the region of the maximum of g. 

P=4 M=3 Spin Glass 
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FIG. 2. Binder parameter as a function of e at tempera- 
ture T = 3.2 for different sizes (L = 3 empty squares, L — 4 
filled circles, L = 5 empty circles, L = 6 filled squares). 
The same plot for the Ed wards- Anderson model (inset) with 
T = 1.3, 1.4, 1.5 {To ^ 1.2 O) shows a completely different 
behavior compared to the p — 4 model. 

In what follows we try to estimate the shape of this 
critical line using finite-size scaling techniques. A de- 
tailed investigation of the critical line as well as its critical 
exponents is presently out of reach due to the smallncss 
of the sizes studied. Still we can approximately deter- 
mine its shape. Let us suppose (as data of figure || sug- 
gests) that there is a critical line e ~ C(T — T C ) A where 
C, A > and T c ~ 2.60. Assuming the validity of the 
scaling hypothesis we may write g(T, e) = g(e(T— T C )~ A ). 
In figure [| we plot the scaling behavior within the scal- 
ing region for different values of T and e for the largest 
size L = 6. The scaling is quite good and proves two re- 



sults: 1) The position of the maximum stays along a well 
defined line in the (e, T) plane and 2) The value of the 
maximum of g is the same everywhere along that line. A 
good collapse of data is obtained with an exponent A ~ 2. 
The position of the maximum of the scaled data yields 
a value of C ~ 0.17 and g ma x — 0.2. This value is uni- 
versal along the critical line and approximately coincides 
with the value of the Binder parameter at the tricritical 
point ||. Although we expect that the critical line will 
have some L dependence, very similar results obtained 
for L = 5 indicate that our estimate of the critical line 
for L — 6 is a good estimate of the true (infinite volume) 
line. 
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FIG. 3. Binder parameter plotted as a function of 
e(T - 2.6)~ 2 in the range .01 < e < .09, 3.0 < T < 4.0 for 
L = 6. In the inset we show the estimated critical line (see 
the text) compared to the points for L = 5, 6 where appears 
a maximum of g in the plane (e, T). 

Now we discuss the physical interpretation of this crit- 
ical line. As already said in the introduction, the pres- 
ence of this transition line is consequence of the exis- 
tence of an exponentially large number (with the volume 
exp(aU), a being a coefficient) of metastable states in 
this model. In mean-field theory (AGM or ideal MCT) 
metastable states have an infinite lifetime so there are 
infinitely large barriers which separate them. This is the 
reason why in mean-field theory ergodicity already breaks 
at the mode-coupling (also called dynamical) tempera- 
ture. In short-ranged systems or real glasses metastable 
states decay by nucleation processes so ergodicity is al- 
ways restored. A typical feature of the ideal MCT sin- 
gularity is the characteristic two-step relaxational decay 
in correlation functions, the so-called a and (3 processes. 
In ideal MCT the typical relaxation time associated to 
the a process diverges at and the ergodicity param- 
eter jumps discontinuously at T^. An accurate study of 
correlation functions reveals that the two-step charac- 
teristic relaxation curve is absent in the present model 
(@) HI1- The absence of a plateau in these curves indi- 
cate that do not exist two well separated time-scales (at 
and /3), like in generic glass- forming liquids, but a con- 
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tinuous hierarchy of time scales for nucleation processes. 
We interpret this result as consequence of the continuous 
nature of the transition everywhere in the critical line 
e — T . A possible scenario for the potential function for 
this type of transition is depicted in figure ^. There is 
not a typical time scale for nucleation process (where a 
small excited droplet or bubble decays from q — qEA to 
q = 0) and the potential around the secondary minimum 
q = Qea is marginally stable. Across the critical line 
the Edwards- Anderson parameter qEA is continuous but 
dq ^ A is discontinuous, the potential being completely flat 
in q. Note that at T = T c , e = where the complexity 
vanishes, the Edwards- Anderson parameter is also con- 
tinuous in agreement with the absence of the two-step 
relaxation in this model. The complexity in this model is 
defined by the height of the secondary saddle point which 
moves towards q = for T = T c . The present scenario 
is very different to that found in the Edwards- Anderson 
model in finite dimensions where no additional spinodal 
point (apart from the paramagnetic one q = 0) is found 
above T c . Concerning dynamical processes and thermal- 
ization effects this model behaves also quite differently. 
Spatial regions may be frozen if their local temperature 
(measured always respect to the intensity of the local 
interaction) is low enough for nucleation processes to de- 
cay very slowly. This is apparent from figure |] where the 
secondary saddle point may become a stable minima in 
certain Griffiths regions. This effect has been observed 
in numerical studies of the P(q) where a secondary peak 
in that distribution function has been observed already 
for small sizes for certain disorder realizations B . 
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model which in the mean-field approximation is well con- 
sistent with the Adam-Gibbs-DiMarzio theory and with 
the ideal mode-coupling theory. In the e — T plane we 
have found evidence for a line of critical points. The 
sole existence of this transition suggests that a bunch 
of metastable states do change the high-temperature be- 
havior through their configurational entropy or complex- 
ity The interest of this study is that the full analysis 
is done in the high-temperature phase where thermaliza- 
tion is easier to achieve (but still, huge thermalization 
times are needed for large sizes) compared to the low- 
temperature region. To our knowledge, such an analy- 
sis of this paramagnetic-paramagnetic transition has not 
been previously done for spin glasses. The main caveat of 
our results is that we cannot exclude the possibility that 
we are observing a crossover behavior where the max- 
imum value of g would eventually vanish for L —> oo 
|l4| . Unfortunately, metastability slows very much ther- 
malization even in the high-temperature region. Further 
studies should extend the present analysis to the study 
of the potential energy landscape and a numerical esti- 
mate of the configurational entropy for this model (such 
as has been done for glass forming liquids |15|) as well 
as the role of the Griffiths singularities in the dynamics. 
Finally because a finite e breaks the q — > — q symmetry 
the existence of this critical line cannot be simply ex- 
plained in terms of a low-temperature spin-glass phase 
characterized by only two symmetrically related equilib- 
rium states. So probably one must invoke more complex 
descriptions for the ground state structure and the type 
of excitations in the spin-glass phase. 
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FIG. 4. Potential function for the short-range p-spin spin 
glass at e — above T c as a function of temperature (from 
high to low temperatures from above to below). The sec- 
ondary minimum at q ^ is a spinodal point. The height of 
that secondary saddle point is the complexity (the logarithm 
of the number of metastable states per site) which vanishes at 
T c . The potential at T = T c becomes flat and qea vanishes 
at T c (following the dashed line). 

In summary, we have studied a short-range spin-glass 
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